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ABSTRACT 

We present new thermal equilibrium solutions for optically thin and optically thick disks 
incorporating magnetic fields. The purpose of this paper is to explain the bright hard state and 
the bright/slow transition observed in the rising phases of outbursts in black hole candidates. 
On the basis of the results of three-dimensional MHD simulations, we assume that magnetic 
fields inside the disk are turbulent and dominated by the azimuthal component and that the 
azimuthally averaged Maxwell stress is proportional to the total (gas, radiation, and magnetic) 
pressure. We prescribe the magnetic flux advection rate to determine the azimuthal magnetic 
flux at a given radius. Local thermal equilibrium solutions are obtained by equating the heating, 
radiative cooling, and heat advection terms. 

We find magnetically supported {(3 — {p gas +frad) /Pm&g < 1), thermally stable solutions for 
both optically thin disks and optically thick disks, in which the heating enhanced by the strong 
magnetic field balances the radiative cooling. The temperature in a \ow-/3 disk (T ~ 10 7 ~ n K) 
is lower than that in an advection dominated accretion flow (or radiatively inefficient accretion 
flow) but higher than that in a standard disk. We also study the radial dependence of the thermal 
equilibrium solutions. 

The optically thin, low-/? branch extends to M > O.lMEdd where M is the mass accretion 
rate and A/Edd is the Eddington mass accretion rate, in which the temperature anti-correlates 
with the mass accretion rate. Thus optically thin low-/? disks can explain the bright hard state. 
Optically thick, low-/3 disks have the radial dependence of the effective temperature T c g cx tu~ 3/ " 4 . 
Such disks will be observed as staying in a high/soft state. Furthermore, limit cycle oscillations 
between an optically thick low-/? disk and a slim disk will occur because the optically thick low-/? 
branch intersects with the radiation pressure dominated standard disk branch. These limit cycle 
oscillations will show a smaller luminosity variation than that between a standard disk and a 
slim disk. 

Subject headings: accretion, accretion disks — black hole physics — magnetic fields — X-ray: binaries 
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1. Introduction 

X-ray spectral states of galactic black hole can- 
didates (BHCs) have been classified on the basis 
of their spectral shape and X-ray luminosity. In 
the high/soft state, the X-ray spectrum is dom- 
inated by a modified blackbody (thermal) com- 
ponent with characteristic temperature ~ 1 keV. 
The low/hard state is observed at low luminos- 
ity (L < O.lLEdd), where LEdd = 4:TtcGM/k cs ~ 
1.47 x 10 39 (M/10M q )(k os /0.34 cn^g" 1 )- 1 erg s" 1 
is the Eddington luminosity. Here M is the black 
hole mass and K es is the electron scattering opac- 
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ity. The X-ray spectrum in the low/hard state 
is dominated by a power-law component (pho- 
ton index T ~ 1.7) with an exponential cutoff 
at ~ 200 keV. The very high/steep power law 
state is observed at high luminosity (~ ^Edd), in 
which the blackbody component and the power- 
law component become comparable. Super crit- 
ical accretion flows have been identified by some 
auth ors in narrow-line Seyfert 1 galaxies (NLSls; 
e.g., iKawaeuchil [2003) and in ultraluminous X- 
ray s ources (ULXs; e.g., IVierdavanti et al. 2006. 
2008). The spectral fitting models used in their 



paper are based on th e slim disk model (e.g., 



Abramowicz et al. 19881 ). These objects are re- 



ferred to as staying in the slim disk state. 

Recently, iMivakawa et al. ( 20081 ) analyzed the 
results of RXTE observations of the black hole 
candidate GX 339-4 in the rising phases of tran- 
sient outbursts. When the luminosity is lower 
than ~ 0.07LEddj the X-ray spectrum is domi- 
nated by a power-law component with an expo- 
nential cutoff, and the cutoff energy is roughly 
constant ~ 200keV. This is a characteristic fea- 
ture of the low/hard state. By contrast, when the 
luminosity becomes higher than ~ 0.07LEdd, the 
cutoff energy anti-correlates with the luminosity 
and decreases to ~ 50keV. This suggests that the 
electron temperature of a disk decreases as the lu- 
minosity increases. Furthermor e, this state has 
been o bserved up to ~ 0.3LEdd- IMivakawa et al 



(2008) have named it the "bright hard state". 

Gierlihski fe Newtonl (|2006l ) have reported two 
types of hard-to-soft transitions. One is the 
bright/slow transition that occurs at ~ 0.3LEdd 
and takes more than 30 days. The other is the 
dark/fast transition which occurs at less than 
0.1.Lehh and takes less than 15 days (see also 
iDone fc Gierlihskill2003l ). 



In the conventional theory of accretion disks, 
the concept of phenomenological a-viscosity is in- 
troduced. In this framework, the tutp-component 
of the stress tensor is assumed to be propor- 
tional to the sum of the gas pressure p gas and 
radiation pressure p rad , t w<p = ass(p gas + Prad)- 
Her e cers is the viscosi t y par ameter introduced 
by IShakura fc Sunvaevl ( 1973). The standa rd 



disk m odel was developed bv lShakura fc Sunvaev 
(|1973l ). under the assumption that a substantial 
fraction of the viscously dissipated energy in an 
optically thick disk is radiated locally. This stan- 



dard disk model reproduced the thermal compo- 
nent in the high/soft state but cannot explain the 
hard X-ray spectrum observed in the low/hard 
state. 



Thorne fc Pried (|1975f ) proposed that hard X- 



rays from Cyg X-l are produced in an in n er op 



tically thin hot disk. IShibazaki fc Hoshil (|l975l ) 



studied the structur e and stability of optically thin 
hot ac cretio n disks. Eardlev. Lightman. fc Shapiro! 
(| 19751) and lShapiro. Lightman. fc Eardlevl (|l976l 
hereafter SLE) constructed a model for optically 
thin two-temperature accretion disks in which the 
ion temperature is higher than the electron tem- 
perature. However, these disks are thermally un- 
stable. 



Ichimarul (| 19771 ) pointed out the importance 
of energy advection in hot, magnetized accretion 
flows, and obtained steady solutions of optically 
thin disks. In this model, the viscously dissipated 
energy is stored in the gas as entropy and ad- 
vected into the central object rather than being 
radiated. Such geometrically thick, optically thin, 
advection-dominated accretion flows (ADAFs) 
or radiatively inefficient accretion flows (RIAFs) 
have be en st udied ex tensively by Narayan fc Yil 



(11994 Il995l). a nd lAbramowicz et all (|l995fT 



Esin et al.l ( 19971 ) found that the maximum lu 



minosity of the ADAF/RIAF solutions is 
0.4a 2 LEdd- These solutions explain the low/hard 
state. 

However, the ADAF/RIAF solutions cannot 
explain the high luminosity observed in the bright 
hard state and during the bright/slow transition 
unless a ~ 1. Furthermore, the ADAF/RIAF 
solutions do not show an anti-correlation be- 
tween the luminosity and the electron tempera- 
ture (therefore, the cutoff energy) observed in the 
bright hard state. 

In the light of this discussion, our main aim 
is to propose steady models of magnetically sup- 
ported black hole accretion disks that can explain 
the bright hard state and the bright/slow transi- 
tion. 

(jl989h 



Okada et al 



constructed a model of 
non-relativistic equilibrium tori with purely az- 
imuthal m agnetic fields in the pseud o-Newtonian 



potential (jPaczvhskv fc Wiita 



1980( ). which has 



been used as in itial conditions of numerical sim- 
ulations (e.g., Machida. Havashi. fc Matsumotol 
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2000; Ma chida fe Matsumotoll2003l ; IMachida et al 



2006). Komissarovl ( 2006 ) developed a model of 
relativistic equilibrium tori with strong toroidal 
magnetic fields around a rotating black hole. How- 
ever, the angular momentum transport was not 
taken into account in their models. 

The origin of the efficient angular momen- 
tum transport (that is, q-visc osity) in accretion 
disks had been a puzzle until Balbus fc Hawlevl 
(|199ll ) pointed out the importance of the magneto- 
rotational instability (MRI). The MRI can excite 
and maintain magnetic turbulence. The Maxwell 
stress generated by turbulent magnetic fields dom- 
inates the total stress and efficiently transports the 
angular momentum of the disk material. Many re- 
searchers have studied the growth and saturation 
of the MRI by conducting local and global mag- 
netohydrodynamic (MHD) simulations. 

Global three-dimensional MHD simulations of 
optically thin, radiatively inefficient accretion 
disks indicate that the amplification of magnetic 
fields becomes saturated when j3 = p gas /Pmng ~ 
10, and that then the disk approaches a quasi- 
steady state. In this quasi-steady state, the ratio 
of the azimuthally averaged Maxwell stress to the 
azimuthally averaged gas pressure, which corre- 
sponds to the a-parameter in the conventional 
theory o f accr e tion disks, is ~ 0.01 — 0.1 (e.g., 
Hawlevl l2000t IMachida. Havashi. fc Matsumotol 
2000HHawlev fc Krolilj200 it IMachida fc Matsumoto 



2003) 



Three-dimensional MHD simulations have re- 
vealed that magnetic fields contribute to not only 
the angular momentum transport but also to 
disk heating via the dissip ation of magnetic en- 



ergy inside the disk (e.g., Machid a et al. 2006; 



Hirose et al.1l2006t iKrolik et alj |2007) 



Machida et al 



(|2006h found that when the 
mass accretion rate exceeds the threshold for the 
onset of a cooling instability, an optically thin, ra- 
diatively inefficient, hot, gas pressure dominated 
accretion disk (ADAF/RIAF-like disk) evolves 
towards an optically thin, radiatively efficient, 
cool, magnetic pressure dominated accretion disk 
(low-/3 disk). In their simulation, when the den- 
sity of the accretion disk exceeds a critical den- 
sity, the cooling instability taking place in the 
ADAF/RIAF-like disk results in a decrease in the 
gas pressure. Thus the disk shrinks in the verti- 
cal direction. During that time, the azimuthal 



magnetic flux (B V )H is almost conserved be- 
cause the cooling time scale is faster than that 
of the buoyant escape of the magnetic flux, where 
(B v ) is the mean azimuthal magnetic field, and 
H is the half thickness of the disk. As the disk 
shrinks, the magnetic pressure exceeds the gas 
pressure because the magnetic pressure increases 
due to flux conservation. Then the disk is sup- 
ported by magnetic pressure and stops shrink- 
ing. In such low-/3 disks, the magnetic fl ux can- 
not e scape through the Parker instability (| Parker 



19661 ) because the strong magnetic tension sup- 



presses the growth of the Parker ins tability (see 
Shibata. Tajima. fc Matsumotol [l990h . Therefore 



the disk stays in a quasi-equilibrium cool state. 

In this state, the dissipated magnetic energy is 
mainly radiated (not advected), and the Maxwell 
stress is proportional to the sum of the gas pres- 
sure and the magnetic pressure. The luminosity 
of the disk becomes higher than that of the initial 
ADAF/RIAF state. These numerical results indi- 
cate that a luminous hard state exists in BHCs. 

lOda et all ( 2007 ) constructed a steady model 
of optically thin accretion disks on the basis of 
these results from three-dimensional MHD simu- 
lations of accretion disks. They found that when 
the mass accretion rate exceeds the threshold for 
the onset of thermal instability, the thermally and 
viscously stable low-/3 branch app ears. These so- 
lution s confirmed the results by IMachida et al. 
(<2006l) . On this branch, the luminosity emitted 



by an accretion disk can exceed 0.2LEdd- Further- 
more, the mean temperature is lower than that of 
the ADAF/RIAF (higher than that of the stan- 
dard disk) and anti-correlates with the mass ac- 
cretion rate. Therefore, they concluded that opti- 
cally thin low-/3 disks can explain the bright/hard 
state. 

In this paper, we extend the model of optically 
thin low-/3 disks to the optically thick regime. We 
obtain local thermal equilibrium curves of accre- 
tion disks incorporating the magnetic fields. The 
total pressure is the sum of the gas pressure, the 
radiation pressure, and the magnetic pressure. We 
consider the thermal bremsstrahlung radiation in 
the optically thin limit and blackbody radiation 
in the optically thick limit as radiative cooling. 
To complete the set of basic equations, we spec- 
ify the advection rate of the azimuthal magnetic 
flux. The basic equations are presented in § [2] In 
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§ [3] we obtain the thermal equilibrium curves and 
study their radial dependence. § [4] is devoted to a 
discussion. We summarize the paper in § [5] 

2. Models and Assumptions 
2.1. Basic Equations 

We extended the basic equations for one- 
dimensional stead y black hole accretion flows (e.g., 
Kato et all 120081 ). incorporating magnetic fields. 
We adopt cylindrical coordinates (vj, (f, z). Gen- 
eral relativistic effects are simulated using the 
pseudo-Newtonian potent ial tp = —GM/(r — r s ) 
(|Paczvhskv fc Wiitallf980h . where G is the gravi- 
tational constant, M is the black hole mass (we as- 
sume M — fOM© in this paper), r = (vj 2 + z 2 ) 1 / 2 , 
and r s — 2GM/c 2 is the Schwarzschild radius. 

We start with the resistive MHD equations 



dp 
Of 



+ v • (H = o , 



(i) 



-pVlp - V (p gas + p rad ) + 



jxB 



showed that magnetic fields inside the disk are 
turbulent and dominated by the azimuthal com 



ponent in a quasi-steady state fe.g. jMachida et al 
20061: iHirose et ai1l2006l ). On the basis of results 



of these simulations, we decomposed the mag- 
netic fields into the mean fields B = (0, (B v ),0) 
and fluctuating fields SB = (SB^^SB^^SBz) and 
also decomposed the velocity into the mean veloc- 
ity v = (v^,v v ,v z ) and the fluctuating velocity 
Sv = (<5u ro , Sv v , 5v z ). Here ( ) denotes the az- 
imuthal average. We assume that the fluctuating 
components vanish when azimuthally averaged, 
(Sv) = (SB) = 0, and that the radial and vertical 
components of the magnetic fields are negligible 
compared with that of the azimuthal component, 
\{B V ) +<5B V | > \SB„\, \SB Z \ (see panel (a) in 
figure [J) . 

Let us derive the azimuthally averaged equa- 
tions. We assume that the disk is in a steady 
state and in hydrostatic balance in the vertical 
direction. By azimuthally averaging equations ([1]) 
- (|4]) and ignoring the second order terms of Sv, 
SBzz, and SB Z , we obtain 



(2) 



1A 

vj dw 



d 

(wpv^) + — (pv z ) = 



d 



(5) 



d(pe) 
dt 



+ V • [(pe + p gas + J5 rad ) V] 

- (« ' V) (p gas + Prad) = 9 + ~ (f ,(3) 



PV U 



' dw 



+ pv z 



-P 



dz w 

di> d Ptot (b 2 ) 

> (6) 



dz 



Amw 



dB 

~dt 



V x [ v x B 



4tt 



-V3 



(4) 



where p is the density, v is the velocity, B is the 
magnetic field, j = cV x B/Air is the current den- 
sity, and pe = p g&s / (7 — 1) + 3p ra d is the internal 
energy. Here, 7 = 5/3 is the specific heat ratio and 
Pgas = T^-pT/ p is the gas pressure where 1Z is the 
gas constant, p is the mean molecular weight (we 
assumed to be 0.617), and T is the temperature. 
In the energy equation ^ , q + is the heating rate, 
and q~ is the radiative cooling rate. In the in- 
duction equation ([3]), 77 = c 2 /4wa c is the magnetic 
diffusivity, where o~ c is the electric conductivity. 



dv v dv v pv^v v 

pv m -^- + P v z-^r- 4 

dw dz 



J__d_ 

w 2 dw 



2 (B m B v ) 
An 



d_ 

dz 



w 

[B V B Z ) 

47T 



= _chp_ _ 1 dp tot 
dz p dz 



~ Vz TT (Pgas +Prad) 

dz 



(7) 
(8) 



—— [{pe + p gas + Prad ) V m \ H (pe + p gas + p ra d ) 

ovj vj 

+ ir[(P e + Pgas + Prad) V Z 

dz 



9 , 

V ™~Q^ (Pg as +Prad) 



(9) 



2.1.1. Azimuthally Averaged Equations 

Three-dimensional global MHD and local Radiation- 
MHD simulations of black hole accretion disks 



oz ovj 
-{V x (Sv x SB)} V - {r/V x (V x B)} v , (10) 
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Pgas+Prad+Pmag is the total pressure 



B 2 ) /8tt is the azimuthally averaged 



where p to t 
and />■■ — , 
magnetic pressure. The third and fourth terms 
on the right-hand side of equation (TIT))) represent 
the dynamo term and the magnetic diffusion term 
which we approximate later on the basis of the 
results of the numerical simulations. Note that we 
have considered only the azimuthal component of 
the induction equation (|4|). 



2.1.2. Vertically Integrated, 
Equations 



Azimuthally Aver- 



We assume that the radial velocity the 
specific angular momentum £ — vjv v , and f3 = 
(.Pgas + Prad)/Pmag are independent of z, and we 
assume a polytropic relation p tot — Kp 1+1 / N in 
the vertical direction, where TV is the polytropic 
index (we adopt N = 3 in this paper). Vertical 
integration of equation (|8|) yields 



p(w,z) = po(ro) \ l -jp 



N 



Ptot(w, Z) = Ptoto(ro) ^1 - 

and we assume 

T(w,z) = T {m) (1-^2 



N+l 



(11) 



(12) 



(13) 



where H is the half thickness of the disk given by 



n^H 2 = 2(N+1) 



PtotO 

Po 



(14) 



where ft K o = {GM/w) 1/2 / (w - r s ) is the Keple- 
rian angular velocity. Here the subscript refers 
to quantities in the equatorial plane. Under these 
assumptions, the surface density S and the verti- 
cally integrated total pressure W to t are given by 

i-H 

pdz 

-H 

Poll 

-H \ 

2 Po I N H , (15) 



z 



dz 



H 



Ptotdz 



H 

PtotO [ 1 

-H 

2ptotolN+iH , 



H 2 J 



N+l 



dz 



(16) 



where I N = (2 N JV!)/(2JV+1)!. Using E and W u 



equation (fT4[) can be rewritten as 



fl^H 2 = {2N + 3) 



(17) 



We note that the half thickness of the disk depends 
on not only the gas and radiation pressure but also 
on the magnetic pressure. 

We now integrate the other basic equations in 
the vertical direction and assume Keplerian rota- 
tion (fi = Qko) instead of equation ©. We obtain 



M 



n = n 



KO 



M{£ K0 - 40 = -2irm 2 



H 



H 



(18) 



(19) 



' B ™ Btp) -dz (20) 



-in 



M W e; 



ad 



2ttw 2 



$ = 



t = Q + - Q 



rad ' 



(21) 



/ v^{B ip )dz 

J-H 



K7out rH 



[{V x (Sv x SB)}^ 



H 



H 

- {77V x (V x B)} v ]dwdz + const. (22) 

where M is the mass accretion rate, £ko = ^ 2 ^ko 
is the Keplerian angular momentum and £- m is 
the specific angular momentum swallowed by the 
black hole. W sas and W ra d arc the vertically in- 
tegrated gas and radiation pressures, respectively. 
The left-hand side of the energy equation has been 
simplified by in troducing the ent ropy gradient pa- 
rameter £ (see Kato et all 120081 ). On the right- 
hand side, Q + is the vertically integrated heating 
rate, and Q~ ad is the vertically integrated radia- 
tive cooling rate. In equation (f22|) . $ is the ra- 
dial advection rate of the azimuthal magnetic flux 
(hereafter we call it the magnetic flux advection 
rate). 



■5 



2.2. a-prescription of the Maxwell Stress 
Tensor 



a disk, and set the vertically integrated heating 
rate as follows: 



Global MHD simulations of radiatively ineffi- 
cient, magnetic press ure supported (low-/?) disks 



( Machida et al. 2006f) showed that the ratio of 



the azimuthally averaged Maxwell stress to the 
sum of the azimuthally averaged gas pressure and 
magnetic pressure is nearly constant except in 
the plunging region (a B = -(5 ro 5 v /47r)/(|5 gas + 
Pmag) ~ 0.05 - 0.1). Global MHD simulations 
of radiatively inefficient, gas pressure supported 
(high-/?) disks also implied that «b is nearly con- 
stant because a m = — (B a B (f ,/4Tr)/(p mag ) ~ 0.2 — 
0.5 a nd ~ 10 are nearly co nstant inside the disk 
(e.g., lHawlev fc Krolikll200lh . 

Following the simulation results, we assume 
that the azimuthally averaged w^j-component of 
the Maxwell stress inside a disk is proportional to 
the total (gas, radiation, and magnetic) pressure 



{B^B V ) 

47T 



(23) 



Integrating in the vertical direction, we obtain 



H 



—H 



[B^B^ 

47T 



-dz — —aWt< 



(24) 



This is one of the key assumptions in this paper. 
When the magnetic pressure is high, the stress 
level can be high even though the gas pressure 
and the radiation pressure are low. 

2.3. Heating and Cooling Rates 

2.3.1. The Magnetic Heating Rate 

In the conventional theory, the viscous heat- 
ing was expressed as q^ is = t mip ru (dfl/dw), where 
t mtp is the wy-component of the total stress and 
57 is the angular velocity, respectively. The re- 
sults of three-dimensional MHD simulations indi- 
cate that dissipation due to the thermalization of 
magnetic energy dominates the total dissipative 
heating rate throughout a disk and is expressed as 



(B^jB^/i^w ( dil/d v u) (e.g. 



20061 : iMachida et all 120061 : iKrolik et al 



Hirose et al.. 
2007h . 



Hereafter, we refer to it as the magnetic heating 
rate. 

Following these simulation results, we employ 
magnetic heating as the heating mechanism inside 



Q + = 



H 



H 



(B m B v ) dn 

TU— — 

4-7T dvo 



dz = -aWtotw^— , (25) 
dvo 



where we have used equation (f2"3")l . We note that 
if the magnetic pressure is high, the heating rate 
can also be large even when the gas pressure and 
the radiation pressure are low. 

2.3.2. The Radiative Cooling Rate 

We assume that radiative cooling is provided by 
the thermal bremsstrahlung emission in the op- 
tically thin limit. Thus the vertically integrated 
optically thin cooling rate is given by 



H 



' —H 



6.2 x 10 20 p 2 T 1/2 dz 



= 6 .2xlO M ^±^$ T o 1/a .(26) 
2//v B 

In addition, we assume radiative cooling due to 
blackbody radiation in the optically thick limit. 
The vertically integrated optically thick cooling 
rate is expressed as 



Q 



thick 



16o-InT$ 
3t/2 



(27) 



where a is the Stefan-Boltzmann constant, r = 
Tabs + t cs is the total optical depth, r a b s is the 
absorption optical depth, r os = 0.5re e s£ is the 
electron scattering optical depth, and k cs = 
0.34 cm 2 g _1 is the electron scattering opacity. 

Although we should solve the transfer equation 
in the intermediate c ase, we use the following ap- 



proximate form (e.g., Hubenvl | l990l : Naravan fc Yi 



Il995t lAbramowicz et al.lll996l ). 

16aI N T 4 



where 



Tabs 



Q 



thin 



3r/2 + V3 + T abs - 



6.2 x 10 20 Isn+i/2 S ; 



Qthick 



16(7 



21 



N 



H 



(28) 



TcT 7/2 ■ (29) 



Using these approximations, the vertically inte- 
grated radiation pressure and the total pressure 
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can be expressed as 
Q- h 



ad 



4c I 



H 



N 



3r 
Y 



73 



(30) 



(i + r 1 ) 



rad "I - ^mag 



I N [i 



T E 



4c 7jv 



,(31) 



respectively. 



2.4. Prescription of the Magnetic Flux 
Advection Rate 

We complete the set of basic equations by speci- 
fying the radial distribution of the mean azimuthal 
magnetic fields. If we ignore the dynamo term 
and the magnetic diffusion term and perform the 
integration in the left-hand side of the induction 
equation (|22|) , we obtain 



$ = -VvjB {m)2 



'(JV+l)/2 



H = const. (32) 



where 



n. 



1/2 



Bo{w) = 27T 1 / 2 (^ToJ Y}' 2 ^ 1 ' 2 ^ 1 ' 2 (33) 

is the mean azimuthal magnetic field in the equa- 
torial plane. However, $ is not always conserved 
in the radial direction due to the presence of the 
dynamo term and the magnetic d iffusion term. 
According to Machida et al. ( 2006f ). this equals 



$ oc vj^ 1 in the quasi steady state. Following 
this result, we parametrize the dependence of $ 
on vj by introducing a parameter, £, as follows. 



*(w;C,M) = *out(M) 



-c 



(34) 



where $ ut is the magnetic flux advection rate at 
the outer boundary vj — Tu out . When Q — 0, the 
magnetic flux is conserved in the radial direction, 
and when £ > (or ( < 0), the magnetic flux in- 
creases (or decreases) with a decreasing radius (see 
panel (b) in figure [1} . The azimuthal magnetic 
flux can increase inside a disk when the azimuthal 
flux of opposite polarity buoyantly escapes from 



the d isk (e.g., Nishikori. Machida fc Matsumotol 
20061) 

Equation (|34j) is the second key assumption in 
this paper. Specifying the magnetic flux advec- 
tion rate enables the magnetic pressure to increase 
when the disk temperature and the disk thickness 
decrease. By contrast, if we specified the plasma (3 
at each radius instead of the magnetic flux advec- 
tion rate, the decrease in the temperature results 
in a decrease in magnetic pressure. This is incon- 
sistent with the r esults of three-dimen sional MHD 
simulations (e.g.. Machida et al.1120061 ). 

3. Results 

3.1. Local Thermal Equilibrium Curves 

We solved the above basic equations at vj = 
5r s to obtain local thermal equilibrium curves. 
We also assume £[ n — £k(5t s ) in the <p compo- 
nent of the momentum equation (|20p . Results 
of the global three-dimensional MHD simulations 
of optically thin bla ck hole accretion flows (e.g., 
Machida et al. 20041 ) indicate that £ = 1 (see also 
Oda et al. I I2007IT This positivity of £ means that 
the heat advection becomes a cooling at a fixed 
radius. In this paper, we adopt £ = 1 in both 
the optically thin regime and the optically thick 
regime. 

To determine the flux advection rate at a fixed 
radius, we need to specify $ ut- Hence, we im- 
posed the outer boundary condition that T out = 
TViriai = (lJ,c 2 /3TZ) (ro out /r s ) _1 and (3 out = 10 at 
c^out = 1000r s . Under these conditions, we ob- 
tained $out by solving equations (J20|), (24|), (j3"Tj) . 
and (f34|) . For given M, a, (, an d vj, we obtained 
local thermal equilibrium solutions. 

Figure [2] shows the thermal equilibrium curves 
plotted in the E vs. M/Mem, T , (3, T e g plane 
for a = 0.05, C = 0.6 (thick solid), 0.3 (dashed), 
(dotted), and —1.8 (thin solid), respectively. 
Here M Edd = L Edd /r] c c 2 = 4nGM/ {ri c K ca c) is the 
Eddington mass accretion rate, T} e = 0.1 is the 
energy conversion efficiency, r e g = 0.5^/k C sKh E 
is the effective optical depth, and kb = 6.4 x 
10 22 poTo~ 3 ' 5 cm 2 g _1 is the opacity of free-free 
absorption, respectively. 

When £ has a large negative value (C = —1-8), 
which means that the magnetic field strength at 
vj = 5r s is negligible, the ADAF/RIAF and SLE 



7 



branches appear in the optically thin regime, and 
the standard and slim disk branches appear in 
the optically thick regime as expected. (Strictly 
speaking, although the low-/? branches discussed 
below also appear in much lower temperature re- 
gions, we do not display them because we assumed 
fully ionized plasmas.) 

For large values of £, we obtained new branches 
in both the optically thin and the optically thick 
regime. We call them "low-/3" branches. On the 
low-/3 branches, the dissipated magnetic energy 
is mainly radiated rather than advected. There- 
fore, the disk temperature and the gas pressure are 
lower than in the ADAF/RIAF (but higher than in 
the standard disk). In contrast, the magnetic pres- 
sure is high because the mean azimuthal magnetic 
field increases as H decreases. The low-/? branches 
consist of an optically thin part and an optically 
thick part. We find that the optically thin low- 
(3 branches extend to M > O.lAfEdd correspond- 
ing to L > O.lLEdd when £ > 0.3. We note that 
the critical mass accretion rate for the existence 
of ADAF/RIAF solutions M C , A is ~ 0.1Af E dd and 
that of the gas pressure dominated standard disk 
solutions M Cy s is ~ 0.05AfEdd- The slim disk solu- 
tion appears when M > 5AfEdd- 

It is worth noting that M oc £ and To oc £~ 2 
in the optically thin low-/3 branches, and that 
M oc £ and T oc E 1 / 2 in the optically thick 
low-/? branches. These relations are derived as 
follows. First let us derive the dependence of 
$out(Af) on M. In the range of M < M E dd in 
which the low-/3 branches appear, W gas 3> W ra d 
at n7 out . From equations (fTT)) . (fT8)) . (120]) . and 
dH]), we find that Wtot, out oc M, S out oc M, 
-Hout ~ constant, and i> ro! out ~ constant. There- 
fore equation ([33)) yields B (w ont ) oc Af 1 / 2 and 
equation (J32J) yields $ out oc Af 1 / 2 . Next, we 
derive the dependence of M on E. When the 
magnetic pressure dominates the total pressure, 
W tot ~ B 2 H oc $ 2 /( Uro 2 iJ) oc M/{v^ 2 H). Since 
equations OH), ([20]), and (JT7[) yield w ro oc M/S, 
W tot oc M, H oc (W to t /£) 1/2 oc (M/S) 1 / 2 , we 
find that M oc £ and H is constant. Finally we 
derive the dependence of To on £ by noting that 
magnetic heating (Q + oc Wtot oc M oc £) balances 
radiative cooling on the low-/3 branches. We find 
from Q~ oc £ 2 T 1/2 that T oc IT 2 in the optically 
thin regime, whereas we find from Q~ oc T 4 /E 
that Tq oc S 1 / 2 in the optically thick regime. 



Figure [3] depicts thermal equilibrium curves for 
the same parameters as in figure [2] but plotted 
in the T vs. M/M^dd plane and the r c g vs. 
M/M Edd plane. On the ADAF/RIAF branch, the 
disk temperature is independent of the mass ac- 
cretion rate and is nearly constant at To ~ 10 11 K. 
In contrast, on the optically thin low-/3 branches, 
the disk temperature anti-correlates with the mass 
accretion rate, M oc T ~ 1/2 , and T - 10 7 ~ n K. 
However, it is not possible to compare these re- 
sults directly with the observational data (e.g., 
the cutoff energy or the electron temperature in 
Miyakawa et al. 20081 ) because the disk tempera- 
ture does not represent the electron temperature. 
We will discuss this issue in § 14.21 

Figure [4] shows the same thermal equilibrium 
curves as in figure [2] but for a = 0.01. Compared 
to the results for a = 0.05, the ADAF/RIAF so- 
lution disappears at the lower mass accretion rate 
(M C>A ~ 0.003M Edd ). 

3.2. The Radial Structure of the Disks at 
Thermal Equilibrium 

In the previous subsection, we obtained local 
thermal equilibrium curves by solving the basic 
equations for various M and fixed parameters, a, 
Q, and w. In this subsection, we study the depen- 
dence of the solutions on vj by fixing the param- 
eters, a, and M. Since we focus on the hard- 
to-soft transition of the BHCs, we assume that 
a disk is in the ADAF/RIAF state correspond- 
ing to the low/hard state when the mass accretion 
rate is low. Thus we choose the ADAF/RIAF so- 
lution when three solutions (ADAF/RIAF, SLE, 
and low-/3 or standard disk solution) are found for 
the same mass accretion rate. 

Figure [5] shows the radial structure of a 
disk for a = 0.05, < = 0.6, M/M Edd = 
8.28 x 10" 3 (thick solid), 0.05(dashed), 0.10(dot- 
dashed), 0.13(dotted), 0.45(triple-dot-dashed), 
and 12.0(thin solid), respectively. 

When the mass accretion rate is low (M < 
O.OlAfEdd), the dissipated energy is advected in- 
ward because radiative cooling is inefficient due to 
the low surface density over the whole disk. The 
whole disk is hot and the gas pressure is dominant, 
that is, the disk is in the ADAF/RIAF state. 

As the mass accretion rate increases (O.OlMEdd ^ 
M < 0.2MEdd), the disk undergoes a transition to 
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the optically thin low-/? disk from the outer radii 
because the mass accretion rate exceeds M c ,a at 
the transition radius. Now let us derive the ra- 
dial dependence of M C ,A- We find from the local 
thermal equilibrium curves that the ADAF /RIAF 
branch in which the heating rate balances the 
heat advection rate (Q + ~ Qadv) intersects with 
the SLE branch in which the heating rate bal- 
ances the optically thin radiative cooling rate 
(Q + ~ Qtiiin) at Mc,A- The gas pressure is domi- 
nant in both branches (W to t ~ W gas ). Using these 
relations, we find from the basic equations that 
M Ci a oc a 2 m 5 fl^l \d£lKo / dm\ 3 (£ko — @in)- Since 
M c .a increases inward in the outer region, the 
transitions take place from the outer radius. We 
note that M C) a has its maximum value around 
10r s for l- m = £Ko(3r s ). Therefore, the innermost 
region (w < 10r s ) undergoes a transition to the 
optically thin low-/? disk at a lower mass accre- 
tion rate than in the slightly further out region 
(10r s < m < 30r s ). If we consider transonic solu- 
tions for which l- m < £koj such a feature may not 
appear. 

When 0.2M E dd < M < 5M Ed d, the whole 
disk becomes an optically thin low-/? disk in 
which the dissipated energy is radiated by the 
bremsstrahlung emission. We note that when 
C = 0.6, an optically thick low- /? disk does not 
appear because the large magnetic flux prevents 
the disk from shrinking in the vertical direction. 
Let us derive the radial dependence of the sur- 
face density and the temperature in optically thin 
low-/? disks. In an optically thin low-/? disk, the 
magnetic pressure is dominant (Wtot ~ W mag ) and 
the magnetic heating rate balances the optically 
thin radiative cooling rate (Q + ~ Q^hi n )- Using 
these relations, we find from the basic equations 
that E oc MoT 1 ^/ 5 - 2 (£ K o - 4i) 3/5 ^ko /5 and 
T oc M- 2 a A zD 6 - 4C (dn K0 /dvj) 2 . In the outer re- 
gion where f^Ko °^ nj -3 / 2 and (£ko — ^in) ~ ^ko oc 
tu 1 / 2 , we find that E cx Mor 1 ^/ 5-11 / 10 and 
T cx M-^o^w 1 -^. 

When the mass accretion rate is high [M > 
5AfEdd), the disk undergoes a transition to a slim 
disk in the inner region because the radial heat 
advection becomes efficient again. We note that 
even in a slim disk the magnetic pressure is com- 
parable to or larger than the radiation pressure 
(W mag > IU rad » W gas ) for C = 0.6. 



Figure [6] shows the results for £ = 0. As the 
mass accretion rate increases (0.01M E dd ^5 M < 
0.2AfEdd), the disk undergoes a transition to a 
low-/? disk from the outer radii, and the whole 
disk becomes a low-/? disk when 0.2A/Edd < M < 
5AfEdd- The low-/? disks consist of two parts, an 
inner optically thick part and an outer optically 
thin part. In an optically thick low-/? disk, since 
Q~ - Qthick K 16fT/ArT 4 / (3r/2) oc T cff 4 , we ob- 
tain the radial dependence of the effective temper- 
ature, T c ff oc tu -3 / 4 . This dependence is the same 
as that in the standard disk. In addition, when 
the mass accretion rate is high (M > 5MEdd), 
slim disks consist of an inner radiation pressure 
dominated region and an outer magnetic pressure 
dominated region. 

An example of an extremely weak magnetic 
field model (£ = —1.8) is shown in figure[7] In this 
model, the inner region of a disk undergoes a tran- 
sition to a radiation pressure dominated standard 
disk while the outer region undergoes a transition 
to an optically thick low-/? disk. As a result, the 
inner region of a disk always stays in the high /? 
state. 

When a = 0.05, the gas pressure dominated 
standard disk solutions do not appear because 
M C ,A, being proportional to a 2 , is higher than 
M c ,s f° r 01 = 0.05. The transition to the gas 
pressure dominated standard disk takes place 
when a is sufficiently small so that M C ,A < A^c.s- 
In figure [3 we show the results for a = 0.01, 
C = -1.8, M/M Ed d = 1 x 10~ 3 (thick solid), 
2.02 x 10~ 3 (dashed), 4.09 x 10~ 3 (dot-dashed), 
2.68 x 10~ 2 (dotted), 1.1 5 (triple-dot-dashed), and 
12.0(thin solid), respectively. As the mass ac- 
cretion rate increases (1 x 10~ 3 AfEdd ^ A/ < 
0.03AiEdd), the disk undergoes a transition to a 
gas pressure dominated standard disk in the inner 
region while an optically thick low-/? disk exists in 
the outer region. 

4. Discussion 

4.1. New Thermal Equilibrium Solution 
connecting an Optically Thin Regime 
and an Optically Thick Regime 

We obtained thermal equilibrium curves for 
optically thin and optically thick accretion disks 
incorporating the mean azimuthal magnetic fields. 
We prescribed that the azimuthally averaged 
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Maxwell stress is proportional to the total (gas, 
radiation, and magnetic) pressure. Consequently, 
the heating rate can be large if the magnetic pres- 
sure is high even when the gas pressure and the 
radiation pressure are low. To complete the set 
of basic equations, we specified the magnetic flux 
advection rate rather than the plasma (3 at each 
radius. Hence, a decrease in temperature via the 
cooling instability results in an increase in mag- 
netic pressure due to the conservation of the az- 
imuthal magnetic flux (B V )H. 

Under these assumptions, we obtained low-/? 
solutions in both the optically thin regime and 
the optically thick regime in addition to the 
ADAF/RIAF, SLE, standard, and slim disk so- 
lutions. In the low-/? disks, the magnetic heating 
balances the radiative cooling. The slopes of the 
thermal equilibrium curves plotted in the X vs. T 
plane in figure [Hand figured] indicate that the low- 
(3 disk is thermally stable. The borders between 
the optically thin and thick parts of the low-/3 
branches depend on £ which specifies the mag- 
netic flux at each radius. As the magnetic flux 
advection rate increases, the magnetic pressure 
becomes higher, therefore the magnetic heating 
rate increases. Hence the equilibrium temperature 
becomes higher, and the surface density decreases. 
Thus, the effective optical depth decreases for the 
same mass accretion rate. This means that the 
optically thin part of the low-/? branch extends to 
a higher mass accretion rate as £ increases. 

We also studied the radial dependence of the 
solutions. The ADAF/RIAF undergoes a tran- 
sition to the low-/? disk from the outer radii as 
the mass accretion rate increases. This is because 
M Cj a at outer radii is lower than its value at inner 
radii. Figure [9] shows schematic pictures of con- 
figurations of accretion disks for various M and 

C- 

Let us discuss the stability of l ow-/? disks. 
Shibata. Taiima. fc Matsumotol ( 1990l ) carried out 



MHD simulations of the buoyant escape of the 
magnetic flux due to the Parker instability. They 
showed that once a disk is dominated by the mag- 
netic pressure, it can stay in the low-/3 state be- 
cause the strong magnetic tensio n suppresses the 



growt h of the Parker instability. iMineshige et al 



from the surface layer will reduce the total az- 
imuthal magnetic flux of the disk. This loss of 
magnetic flux can balance the amplification of 
magnetic fields due to the MRI in this layer. The 
results of three-dimens i onal g lobal MHD simula- 
tions bv lMachida et al. ( 20061 ) indicate that low-/? 
disks stay in a quasi-equilibrium state at least for 
a thermal time scale. The longer time scale evo- 
lution of the low-/3 disk is a subject to be studied 
in near future. 

4.2. Bright Hard State undergoing Bright/Slow 
Transition 

One of the main aims of this paper is to explain 
the bright hard state observed in rising phases 



of th e outbursts of BHCs (e.g., iMivakawa et al 



2008). Figure [3] shows that the optically thin part 



(|1995l ) suggested that such low-/3 disks emit hard 
X-rays. When f3 ~ 1 in the surface of the low- 
(3 disk, the buoyant escape of the magnetic flux 



of the low-/? branch extends to > 0.1A/Edd for 
C > 0.3. In the ADAF/RIAF branches, the disk 
temperature is independent of the mass accretion 
rate and remains nearly at the virial temperature. 
This is consistent with the fact that the cutoff 
energy observed in the low/hard state of BHCs 
is independent of their luminosity. In contrast, 
the disk temperature roughly anti-correlates with 
the mass accretion rate in the optically thin low-/3 
branches. This indicates that the electron tem- 
perature in the hard X-ray emitting optically thin 
low-/? disks anti-correlates with the luminosity. 
This feature is consistent with the anti-correlation 
observed in the bright hard state of BHCs. 

In addition, the optically thin low-/? disks 
explain the bright /s l ow tr ansition reported by 
iGierlihski & Newtonl (j2006f l For small C, the 
ADAF/RIAF directly undergoes a transition to 
the optically thick disks emitting the soft X-rays 
at - 0.1M E dd for a = 0.05 (panel (h) -> (g) -> 
(f) in figure El) . This corresponds to the d ark/fast 
transition (see lGierliriski fc Newto"nll2006h . In con- 
trast, for large (, the ADAF/RIAF first undergoes 
a transition to the optically thin low-/? disk, which 
can be identified as the bright hard state. Sub- 
sequently, this disk undergoes a transition to the 
optically thick disk at a mass accretion rate higher 
than the rate for small (. This transition corre- 
sponds to the bright/slow transition, (panel (h) 
-» 0) -» (i) -» (e) in figure [9]) 

It is anticipated that if the magnetic flux es- 
capes from the disk due to magnetic buoyancy or 
dissipates due to magnetic reconnection, the equi- 
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librium curve may shift to the one for smaller £. In 
such a case, the optically thin low-/? disk immedi- 
ately undergoes a transition to the standard disk 
or the optically thick low-/? disk. This transition 
may accompany the ejection of jets observed dur- 
ing the hard-to-soft transition of BHCs because 
magnetic energy stored in the low-/? disk is re- 
leased. 

Although we have adopted the single tem- 
perature model (Tj = T ), it has been pointed 
out that the electron temperature becomes lower 
than the ion temperature (Tj ~ 10 11_12 K) in 
such a low density, high tem perature region (e.g., 



low-/? disk (panel (g) in figure [9]). An increase in 
seed photons can result in a decrease in the tem- 
perature of the Compton corona. Furthermore, 
the hard X-rays Comptonized in the corona can 
be incident on the optically thick low-/? disk and 
can be partly reflected. As a result, the existence 
of such region will enhance the Compton continua 
with a reflection component in the X-ray spec- 
trum. 

To sum up, the optically thick low-/? disk not 
only provides seed photons to the Compton corona 
but also produces a soft spectral excess and a re- 
flection comp onent observed in the low /hard state 



Nakamura et al. 19961 . 19971 ) . Moreover, cooling of Cyg X-l dMakishima et al 



via the synchrotron emission and the inverse- 
Compton effect can be efficient. By contrast, the 
electron temperature and the ion temperature co- 
incide (Tj ~ T c ~ 10 6_7 K) in low temperature 
and high density disks such as optically thick 
low-/? disks, standard disks, and slim disks. In 
subsequent papers, we will take into account the 
decoupling of ions and electrons in the low density, 
high temperature region and incorporate the syn- 
chrotron emission and the inverse-Compton effect, 
in addition to the bremsstrahlung emission as the 
cooling mechanism. 

4.3. Magnetically Supported Thermal Disks 

We obtained solutions for the optically thick 
low-/? disks emitting blackbody radiation (panel 
(f) in figure [9]). The radial distribution of the ef- 
fective temperature in the optically thick low-/? 
disk is roughly the same as that for the standard 
disk (T e g oc tu~ 3 / 4 ) because it does not depend 
on whether the gas pressure or the magnetic pres- 
sure is dominant as long as the heating rate bal- 
ances the blackbody cooling. This indicates that 
the optically thick low-/? disk will be observed as 
a thermal disk. 

The seed photons incident on the Compton 
corona (the low density high temperature region, 
that is, the ADAF/RIAF or the optically thin low- 
/? region) have been thought to be provided by a 
cool disk or by synchrotron emissions. The opti- 
cally thick low-/? disk can also be the source of seed 
photons. When ( < 0, the ADAF/RIAF under- 
goes a transition to the optically thick low-/? disk 
from the radially outer region as the mass accre- 
tion rate increases. Thus, the inner ADAF/RIAF 
region is surrounded by the outer optically thick 



J1655-40 ([Takahashi et al.ll2008f) at the Suzaku X 



l2f)f)8h and GRO 



ray observatory. 

4.4. Limit Cycle Oscillation Between an 
Optically Thick Low-/? Disk and a 
Slim Disk 

So far, we have focused on the hard-to-soft tran- 
sition of BHCs. Now we discuss X -ray flares as 
obser v ed in GRS 1915+105 ( e.g. , iBelloni et af 



19971 : iTaam. Chen fc Swankl Il997l : IPaul et al 



1998; IBelloni et all 12000) . This object frequently 



exhibits state transitions and quasi-regular X- 
ray flares. It has been thought that the lat- 
ter effect may be driven by limit cycle oscilla- 
tions between the standard disk and the slim 
disk (panel (a) (b) in figure ©. lOhsuga 



disk [us 
(2006, 2 



20071) reproduced such limit cycle evolu 
tion by global two-dimensional radiation hydro- 
dynamic (RHD) simulations of optically thick ac- 
cretion disks using the flu x-limited diffusion (FLD; 
iLevermore fc Pomraning 1981 ) approximation for 
the radiation transfer and the phenomenological a 
viscosity for the angular momentum transfer and 
the viscous heating. However, such RHD simula- 
tions do not follow the formation of the low-/? disk 
presented in this paper. 

We first consider a transition from a slim disk 
to an optically thick low-/? disk. We suppose that 
the disk is initially in the slim disk state in which 
magnetic heating balances the heat advection (the 
photon trapping) and that these are greater than 
the radiative cooling (Q + ~ Q a dv > Q ) ■ Fig- 
ure [2] shows that as the mass accretion rate de- 
creases, the slim disk branches disappear. At this 
point, the magnetic heating (and the heat ad- 
vection) become less than the radiative cooling 
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(Q + ~ Qadv < Q~)- Therefore cooling instability 
occurs. As a result, the disk shrinks in the vertical 
direction along the cooling timescale and under- 
goes a transition to an optically thick low-/3 disk 
(panel (e) — > (f) in figure[5]). This situation is sim- 
ilar to the transition from the ADAF/RI AF to the 
optica lly thin low-/? disk found by iMachida et al.l 
( 20061 ). Subsequently, as the mass accretion rate 
increases, the disk undergoes a transition back 
to a slim disk because the optically thick low-/? 
branches disappear (panel (f) — > (e) in figure [9]). 

We emphasize that an optically thick low-/? disk 
undergoes a transition to a slim disk at higher 
mass accretion rates (~ A^Edd for a = 0.5 and 
Q = 0) than those at which a standard disk un- 
dergoes a transition to a slim disk (~ O.lMEdd)- 
This results in a small variation in luminosity for 
quasi-regular X-ray flares compared to that ex- 
pected from the limit cycle between a standard 
disk and a slim disk. 

5. Summary 

We have obtained low-/? solutions for optically 
thin and optically thick disks incorporating the 
mean azimuthal magnetic fields. The key assump- 
tions are (1) we prescribed the Maxwell stress to 
be proportional to the total (gas, radiation, and 
magnetic) pressure, and (2) we specified the mag- 
netic flux advection rate instead of the plasma /? at 
each radius. As a result, the strong magnetic fields 
mainly contribute to the heating and increase the 
half thickness of the disk. In the low-/? disks, the 
heating enhanced by strong magnetic pressure bal- 
ances the radiative cooling rather than the heat 
advection. Optically thin low-/? disks explain the 
bright hard state undergoing the bright /slow tran- 
sition observed in BHCs. The optically thick low- 
(3 disk can be observed as a thermal disk, which 
can supply seed photons into the Compton corona 
(the ADAF/RIAF or the optically thin low-/? re- 
gion) and reflect incident photons from the Comp- 
ton corona. We pointed out the possibility of new 
limit cycle oscillations between an optically thick 
low-/? disk and a slim disk, which will exhibit a 
smaller variation in luminosity than that expected 
from a limit cycle between a standard disk and a 
slim disk. 
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Fig. 1. — Schematic picture of magnetic field lines inside the accretion disk, (a) The magnetic fields 
are decomposed into the mean azimuthal magnetic fields and fluctuating fields. The azimuthal average of 
the radial component of the magnetic fields, is small because positive and negative B CT cancel out. 

However, the azimuthal average of the product of the radial and azimuthal component (B^B^) has a large 
negative value because B^B^ does not change sign when magnetic fields are deformed by nonlinear growth 
of the MM. (b) The advection of the azimuthal magnetic flux. The radial dependence of the magnetic flux 
advection rate <E> is parameterized by £. When £ > 0, $ increases with a decreasing radius. However, $ 
decreases when £ < 0. The mean azimuthal magnetic flux is conserved in the radial direction when C, = 0. 
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log S [ g/cm 2 ] 

Fig. 2. — Local thermal equilibrium curves of accretion disks at w = 5r s on the £ vs. M/MEdd, T, [3 
and r e ff plane for a = 0.05, £ = 0.6(thick solid), 0.3(dashcd), O(dotted), and — 1.8(thin solid, extremely 
weak magnetic field case). MEdd = ^Edd/^eC 2 is the Eddington accretion rate, where the energy conversion 
efficiency is taken to be rj c — 0.1. The low-/? branches connect the optically thin and thick regimes. 
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Fig. 3. — Same as figured] but plotted in the T vs. M / M-&dd plane and the t c q vs. M/MEdd plane. The 
temperature T anti-correlates with M on the optically thin \ow-f3 branches. 
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Fig. 4. — Same as figure [21 but for a — 0.01. 
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Fig. 5. — The radial dependence of (a), temperature T, (b) effective optical depth, r e ff, (c) (3 = {pg&s + 
Prad)/.Pmag) (d) surface density E, (e) the ratio of the heating rate and the radiative cooling rate, the ratio of 
the vertically integrated total pressure and (f) the gas pressure W gas , (g) the magnetic pressure W mag , and 
(h) the radiation pressure W rad for a = 0.05, C = 0.6, M/M Edd = 8.28 x 10~ 3 (thick solid), 0.05(dashed), 
O.lO(dot-dashed), 0.13(dotted), 0.45 (triple-dot-dashed), 12 (thin solid). We chose the ADAF/RIAF solution 
when we obtained three solutions (ADAF/RIAF, SLE and low-/? or standard disk solutions) for the same 
M. 
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Fig. 6. — Same as figure 03 but for £ = 0. 
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Fig. 8.— Same as figure but for a = 0.01, C = -1-8, M/Mmd = 1-00 x 10- 3 (thick solid), 2.02 x 
10- 3 (dashed), 4.09 x 10~ 3 (dot-dashed), 2.68 x 10~ 2 (dotted), 0.45(triple-dot-dashed), 12(thin solid). 
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Fig. 9. — Schematic pictures of the configuration of accretion disks for various mass accretion rates and 
magnetic flux advection parameter £. The optically thin region is indicated by gray, and the optically 
thick region is indicated by black. The panels (c), (d), (g), and (h) represent low mass accretion states 
corresponding to the low/hard state. The panels (i) and (j) represent mildly high mass accretion states 
for large £ corresponding to the bright hard state. The panels (b) and (f) represent the mildly high mass 
accretion state for small C corresponding to the high/soft state. The slim disk state is illustrated in the top 
panels (a) and (e). The left column shows the configuration predicted by the conventional (no magnetic) 
model. 
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